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ABSTRACT 

We show how the second sequence seen lying above the main sequence in cluster 
colour magnitude diagrams results from binaries with a large range of mass ratios and 
not just from those with equal masses. We conclude that the presence of a densely- 
populated second sequence, with only sparse filling in between it and the single star 
main sequence, does not necessarily imply that binary mass ratios are close to unity. 

Key words: stars: clusters of - stars: binary 



1 INTRODUCTION 

With the recent surge of high quality colour-magnitude di- 
agrams (hereinafter CMDs) of clusters obtained both from 
space (Rubenstein & Bailyn 1997, Richer et al. 1997 and 
Elson et al. 1998) and ground based (Ferraro et al. 1997) 
observations it is timely to reconsider the effect of binary 
stars on the observed stellar main sequence. It is well known 
that an unresolved binary system comprising two identical 
stars has the same colour but twice the luminosity of an 
equivalent single star and that such a system, comprising 
two equal-mass main-sequence stars, will appear in the clus- 
ter CMD displaced vertically by 0.753 magnitudes irrespec- 
tive of the wavelength bands used (Haffner & Heckmann 
1937). Because of this, the clear second sequences displaced 
by this amount above the main-sequence visible in many 
cluster CMDs are taken to indicate a population of equal- 
mass binaries (Bergbusch et al. 1991, Kaluzny & Rucinski 
1995 and Santiago et al. 1996) . The fact that the region 
between the single-star main sequence and this second se- 
quence is not very densely populated is then often taken 
to mean that the mass-ratio distribution in the binary sys- 
tems is biased towards equal masses. However it turns out 
that this need not be the case. A system with two unequal 
main-sequence components has a combined colour that is 
redder than the colour of the brighter component as well as 
a luminosity greater than the single star but less than the 
corresponding equal-mass binary (Bolte 1991 and Romani & 
Weinberg 1991). Such a system is displaced both upwards 
and to the right relative to the main-sequence position of 
the brighter component. Thus if we consider the position 
of a system with a mass ratio slightly different from unity, 
relative to the equal-mass system it moves to fainter mag- 
nitude but also to redder colour. Now, because the intrinsic 
slope of the main sequence is such that less massive stars are 
fainter and redder, it is possible for systems with unequal 



components to follow the second sequence downwards and 
to the right in the CMD. We show (section 2) that, for cer- 
tain mass ranges and choice of colours, systems with quite 
extreme mass ratios still lie close to the second sequence, 
0.753 mag above the actual main sequence, appropriate to 
their combined colour. We then show (section 3) how an un- 
biased mass ratio distribution can lead to the clearly sepa- 
rated main and second sequences observed in cluster CMDs. 



2 COMBINING THE COLOURS AND 
MAGNITUDES 

For a given colour index X a main-sequence star of bolo- 
metric luminosity L has an absolute magnitude 



M x = M x ,e ~ 2.5 log 10 — + (5x (T eff ) 



(1) 



where j3x is the bolometric correction appropriate to a main- 
sequence star with effective temperature T e ff and Mx,q and 
Z/0 are the absolute magnitude and bolometric luminosity of 
the Sun. Given a second colour Y a corresponding absolute 
magnitude My can be defined in a similar way and a colour 
by 



(X — Y) = M x ~ My 



(2) 



In the case of an unresolved binary the two components con- 
tribute differently to Mx and My because of their different 
effective temperatures and consequently different bolometric 
corrections. Using the suffices 1 and 2 for the components 
and 3 for the system we can write the combined magnitudes 



M x ,3 = M X , Q - 2.5 log 10 
where 



Lx,i + Lx;. 



(3) 
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Figure 1. Theoretical zero-age single star and equal-mass binary 
main sequences for a metallicity of Z = 0.02 covering a range 
of My. Also plotted for a range of primary masses, Mi, are bi- 
nary points with M2 = <jMi where q ranges from 1.0 — * 0.0 in 
increments of 0.1. The point at q = 0.5 is an open square. 



The combined colour is then 
{X-Y) 3 = M X)3 -M Y , 3 . 



(4) 



(5) 



For a binary system of total mass M = Mi + M% and 
mass ratio q — M2/M1, where we choose the primary mass 
Mi > M2 (corresponding to Li > L2 on the main sequence) 
and hence < q < 1, we can calculate both Li and T ff,i by 
means of the zero-age main-sequence formulae constructed 
by Tout et al. (1996). We can then use the synthetic stellar 
spectra computed by Kurucz (1992) to find the appropriate 
bolometric corrections (5x,i needed to generate a theoretical 
CMD. Figure 1 shows such a CMD containing the zero-age 
single star and equal-mass binary sequences for solar metal- 
licity. Also shown for a range of primary masses are the locus 
of binary points with q ranging from TO — > 0.0 in increments 
of 0.1. We see that, as M2 is decreased for a particular Mi, 
the binary becomes fainter but also redder for q greater than 
some 5crit- Thus the locus of binary points follows the equal- 
mass sequence for mass ratios in the range q CI i t < q < 1. For 
Mi > 3M© we find q cr it ^ 0.5 while for Mi < 2M© we 
find q CI it — 0.9. These ideas were originally formulated by 
Haffner & Heckmann (1937) and later developed by Bettis 
(1975) and Dabrowski & Beardsley (1977) but neither of the 
latter two make any connection between magnitude differ- 
ence and mass ratio. Our results, based on up to date stellar 
models and bolometric corrections, represent a thorough re- 
working. 



3 CLUSTER COLOUR MAGNITUDE 
DIAGRAMS 

To understand how these effects manifest themselves in ob- 
served CMDs we consider a distribution of 5000 stars of 
which 4000 are in unresolved binary systems. We select the 
masses of the single stars and of the binary primaries accord- 
ing to the initial mass function derived by Kroupa, Tout & 
Gilmore (1993) by means of the generating function 



Mi 



0.08 



0.19X 1 55 + 0.05X 
(1.0 -X) ' 58 



(6) 



where we choose X uniformly distributed to give masses in 
the range 1 < Mi/M© < 6. We then choose the mass ra- 
tio in the binary systems so that q is uniformly distributed 
between and 1 to find the secondary masses 



M 2 = qMi 



(7) 



For each single star and binary system we calculate the 
colours and magnitudes as described in the previous sec- 
tion and plot them in the CMD shown in Figure 2. We see 
that the densest population is close to the single star main 
sequence even though most stars are binary. This is because 
faint secondaries make no contribution. With this particu- 
lar mass ratio distribution we see a clear second sequence for 
(B — V) < 0.1 corresponding to Mi > 2M Q . This is much 
clearer in the expanded region shown in Figure 2b than in 
the full CMD in Figure 2a and we note that this is due to 
the size of the symbols used to plot the stars relative to the 
space between them and care should be taken when interpre- 
tating observations by eye. This point CMD is in the form in 
which observations have been presented so far but we note 
that a density contour plot would be less open to false in- 
terpretation (see Figure 3). Between 0.1 < (B — V) < 0.3 
no clear second sequence is visible but for (B — V) > 0.3 it 
is again apparent in Figures 2c and 3c. 



4 DISCUSSION 

We have shown that not only equal mass main-sequence bi- 
nary systems lie on the second sequence 0.753 mag above the 
single star sequence but also many with unequal masses. In 
particular for B and V colours and primary masses Mi > 2 
even systems with quite extreme mass ratios, as low as 
q = 0.5, lie on the second sequence. We have also shown how 
a cluster with uncorrelated binary masses can show a dis- 
tinct second sequence in the CMD. However we must com- 
ment on what we mean by uncorrelated. For Figures 2 and 3 
we have selected each primary mass Mi from the IMF and 
then chosen q < 1 uniformly distributed. Note that choos- 
ing q max > q > 1 uniformly would lead to quite a different 
distribution and it can be argued that we really should be 
choosing > log q > log q m i n uniformly. Alternatively Eggle- 
ton, Fitchett & Tout (1989) denned uncorrelated to mean 
the two masses chosen independently from the same initial 
mass function. For a nonuniform initial mass function such 
a choice does not lead to a uniform distribution for q overall 
nor for any particular primary mass. Indeed, because any 
star is most likely to have a low-mass companion, systems 
with massive primaries will tend to have extreme mass ra- 
tios while those with low-mass primaries will tend to have 
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Figure 2. A) CMD for a randomly generated zero-age population of binary and single stars. B) and C) show two enlarged portions 
of the lower panel which correspond to approximate mass ranges 2 < M\ < 4 and 1 < Mi < 1.4 respectively, both of which show 
well-defined second sequences. 



mass ratios closer to unity. In this sense our uniform distri- 
bution in < q < 1 actually corresponds to a correlation 
in the masses for primaries with Mi > 2 Mq . A CMD with 
masses chosen independently from the initial mass function 
given above does not have a clear second sequence around 
3Mq where extreme mass ratios are more likely but does 
at low masses where systems tend to have equal mass com- 
ponents. To distinguish between the possible mass ratio dis- 
tributions we note that we can employ statistical methods 
such as those described by Kroupa and Tout (1992) to in- 
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vestigate the binary content of Praesepe. In that work the 
colours and magnitudes were implicitly combined according 
to the prescription of section 2. 
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Figure 3. Number density plot corresponding to the same regions as the CMD of figure 2. 1000 000 randomly generated binaries are 
considered with no single stars. The grid is 200 X 200 with 10 contour levels evenly spaced between 0.05 and 0.95 after normalisation of 
the data to the maximum density. 
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